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Abstract 

We studied the pair production rate of charged Higgs bosons 
via gluon-gluon collisions in the general two-Higgs-doublet 
model. The analysis of its production are made with possible 
parameter values of Higgs sectors in this model. We find 
that the production cross sections may be 4 to 41 fermi- 
bar in the future LHC pp collider with these parameters. 
The calculation also shows that this pair production process 
is not sensitive to the neutral Higgs masses, therefore its 
measurment could not give more imformations about the 
neutral Higgs bosons. 
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1. Introduction 



In the sense of Minimal Standard Model(MSM) a neutral Higgs boson is predicted 

without giving its mass value. Therefore the search for Higgs boson is undoubtedly one of 
the important issues in testing the electroweak theory. Untill now, the top quark mass is 
given as 0: m t = 176 ± 10l^ GeV. The absence of Higgs-boson signal events leads to a 
lower bound M H o > 63.5GeV at 95% C.L. obtained 0]. In order to investigate the sectors 
of the theory about which we know very little, the future multi-TeV hadron colliders such 
as the CERN Large Hadron Collider (LHC) and Next Linear Collider(NLC) are elaborately 
designed to devote to the precise characters of Higgs sector and top quark. Besides providing 
the first evidence of Higgs mechanism, the multi-TeV colliders offer the possibility of carrying 
on the research of other unknown particles in the hitherto unexplored mass range. 

So far the extension of the MSM which includes a second Higgs doublet is far from being 
rules out experimentally. One strong motivation for extending Higgs sector into two-Higgs- 
doublet model (THDM), comes from that it has a common feature of many extended models 
of the MSM, such as Minimal Supersymmetric Standard Model(MSSM). The MSSM model 
requires two SU(2) doublets to give masses to up and down quarks[|]. It is therefore of 
interest to determine whether THDM or MSSM is compatible with experiments, and what 
further extensions should be if it is not. 

In future hadron colliders, several mechanisms can produce Higgs boson. Among them 
the Higgs pair productions via gluon-gluon collision, quark-antiquark annihilation or others, 
play a significant role in studing Higgs sectors. The on-mass shell charged Higgs boson pair 
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production with their subsequent cascade decays will provide a good opportunity of studying 
the two Higgs doublet model and a copious source of charged Higgs bosons. Of cause, the 
Higgs pair production has less rate comparing with the single production mechanism, but the 
Higgs pair production via gg and qq annihilation involves the trilinear and quartic couplings 
of the physical Higgs bosons. So the measurement of these couplings may be a good test 
of our theory. At the mean time, the pair production of Higgs boson in hadron collisions 
has a better signal-to-background ratio than one Higgs production process. Therefore it is 
worthwhile to study the precise features of the Higgs pair production processes at LHC, 
NLC and other future colliders which are designed operating in 77,e + e~, pp and pp collision 
modes. Especially when charged Higgs boson mass is not very heavy, its pair production will 
be more significant. Recently T. Plehn, M. Spira and P.M. Zerwas provided the calculations 
of pair production of neutral Higgs particles in gluon-gluon collisions in MSSM0 and charged 
Higgs boson pair production at NLC in 77 modes are also investigated in references JT]. 

In this paper we represent the calculation of the charged Higgs boson pair production 
via gluon-gluon collisions at LHC in the frame of general THDM. The paper is organized as 
follows: The analytical evaluations are given in section 2. In section 3 there are numeriacal 
results, discussion and a short summary. Finally the explicit expressions of the relevant form 
factors are collected in Appendix. 

2. Calculation 

The general THDM contains six free parameters. Generally the four mass parameters, 
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i.e. the parameters for two scalar neutral Higgs boson masses m^o and m H o, one neutral 
pseudoscalar Higgs mass m^o , one pair charged Higgs boson masses m H ± and the additional 
mixing angle a and (3, are chosen as free parameters. The parameter a is defined as the 
mixing angle of the scalar Higgs bosons h° and H°. The angle (3 is a mixing parameter of 
charged and pseudoscalar Higgs bosons and is related directly to the vacuum expectation 
values of Higgs doublets (i.e. tan/5 = v^/v-i). In order to demonstrate the numerical features 
of the charged Higgs boson pair production, we specify a = (3 in our numerical calculation 
for simplicity. 

The generic Feynman diagrams contributing to the subprocess gg — > H + H~ which are 
involved in frame of THDM are depicted in figure 1. The relevant Feynman rules can be 
found in references ||. Since we neglect the couplings between light-quarks and Higgs bosons, 
there are only the following three kinds of production mechanisms: One way is by exchanging 
a 7 or a Z° boson to produce a charged Higgs boson pair (see Fig.l(a,b)), we call it as 7 — Z° 
s-channel. The 7 exchange s-channel diagrams have no contribution and only the axis-vector 
part of the coupling between Z° and the heavy quark in Z° exchange diagram, contributes 
to the cross section. The second mechanism is through producing the virtual neutral Higgs 
h° or Hq bosons which is coupled to gluons by the usual heavy-quark triangle, and the 
charged Higgs boson pair is appeared by the subsequent decay of the virtual Higgs bosons 
(see Fig.l(c,d)), which are called the h° — H° s-channel. Another mechanism is named as 
heavy-quark box diagrams(see Fig.l(e ~ j)). The box-diagrams of exchanging the two initial 
gluons are not plotted in Fig.l. All the box diagrams can be divided into two groups. They 
are corresponding to t-channel and u-channel groups respectively. In this work, we used 
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n-dimensional regularization scheme in the one-loop diagram calculation. 

We denote 9 as the scattering angle between the one of the gluons and final H + boson. 
Then in the system of center-of-mass (CMS) we can express all the four-momenta of initial 
and final particles by means of the total energy V§, the scattering angle 6 and the velocity 
of final charged Higgs bosons. In the CMS of the initial two gluons, the four-momenta of 
the final charged Higgs bosons read in term of following forms 

Pi = (-g-j-g-fe sm6,0,—(3 H ± cos#), 

P2 = — ^/3//± sin#,0, -^p-pH±cos6). (1) 
where the velocity of charged Higgs bosons in the CMS, f3 H ± ) can be represented as 

(3 H ± = ^l-Am 2 H± /s. (2) 

Pz and £> 4 are the four-momenta of the initial gluons and expressed as 

P3 = (—,0,0,—), 

P4 = (^-,0,0,-^-). (3) 
The corresponding matrix element for all the diagrams in figure 1 is written as 

M = M s + M box{i) + M hox{il) 
~ l -,T a T b g 2 s t,{p z )e u { Pi ) 



16n 2 



/ne ra/ VaP2,/3P3, 7 p£ + fue^Pi^m,^)- (4) 

The upper indexes s, box(t) and box(w) represent the amplitudes corresponding to the 
s-channel diagrams, including all the 7 — Z° and h° — H° diagrams, t-channal and u-channal 
box diagrams in Fig.l respectively. The g s is the strong coupling constant and T a (a = 1 ~ 8) 
are the SU(3)c generators introduced be Gell-Mann[g]. The variables s, t and u are usual 
Mandelstam variables in the gluon-gluon center-of mass system with the definitions 

s = ( Pl + P2 ) 2 = (p 3 +p 4 ) 2 , 

t= ( pi -p 3 ) 2 = (pjj -p 4 ) 2 , 

M = (Pi - P4) 2 = (P2 - Ps) 2 - (5) 

We collect the explicit expressions of the form factors /j(z = 1 ~ 12) which appear in 
equation (4) in Appendix. The total cross section for subprocess gg — > H + H~ can finally 



be written in the form 



* + • 2 



^ ) = I6^/- (6) 



where |M| is the initial spin and color averaged matrix element squared and t = (m H ± — 
|s) ± \s(3h±- The color factor including in \M\ is equal to 1/32. The total cross section 
for the charged Higgs pair production through gluon-gluon fusion in proton-proton collisions 
can be obtained by integrating the a over the gluon luminosity. 

v(pp^ gg -> H + H- + X) = f 1 dr^^aigg ^ H + H- at s = rs). (7) 

Jim? ,/s dr 
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Here we introduced a variable r = xi%2 and x\ and xi are defined as 

P3 Pi 

where Pi and P 2 are the four momentum of both incident protons respectively, yfs and 
are the proton-proton and gluon-gluon center-of-mass energies respectively, which can be 
expressed as 

S = (Pi + P 2 ) 2 , S = (p 3 + p 4 ) 2 = Sl^S = TS. (8) 

The quantity of the gluon luminosity can be calculated by using the following equation. 



dr 

dL„q f 1 dX\ 



2 / P 9 (xi,/i)P 9 (-,/i) 

J T Xa l X~\ 



dr Jt Xi 



(9) 



X\ 

The numerical resuts of the process pp —>■ gg — > H + H~ are represented by adopting the 
MRS set G parton distribution function F g (x) ||. 

3. Numerical results and discussion 



In the numerical evolution we take the input parameters [|10j as m& = 4.5 GeV^, 



91.1887 GeV, m w = 80.2226 GeV, G F = 1.166392 x 10~ 5 (GeV)- 2 and a = j^km- The 
top quark mass is just fixed to be the central value m t = 175 GeV. We also used a simple 
one-loop formula for the running strong coupling constant a s as 

a 'W = I - 33-2n f , r, 7 \ ' ( 10 ) 

where a s (mz) = 0.117 and n./ is the number of active flavors at scale /1. The ratio of the 
two vacuum expectation values tan /3 is chosen to be varied in our calculation. The coupling 



strength of the charged Higgs boson with quarks is related to scalar and pseudo-scalar parts, 
each part has factor g + or g_(The definitions of g + and g_ can be found in Appendix). 
These two parts of Yukawa coupling yield a crucial contribition alternatively depending on 
variate tan/5. Because the present experimental data disfavor the small values of tan/3 in 
general THDM. Therefore we limit its quantity in the representative range of 1.5 to 30 in 
the calculation. 

The cross sections of the subprocess gg — > H + H~ with m#± = 150 GeV are depicted in 
Fig. 2(1), on the conditions of m h o = 100 GeV and m H o = 150 GeV. Fig. 2(2) shows also the 
cross sections of this subprocess as the functions of CMS energy of incoming gluons with the 
same neutral Higgs boson masses, but m H ± — 300 GeV. In both figures, the full line curves 
and dashed line curves are corresponding to the tan/3 = 1.5 and tan/3 = 30 respectively. 
These graphs show that the cross section of the subprocess is approximately two times larger 
when the ratio of vaccum expectation values is 30 than that when tan/3 = 1.5 in the energy 
region less than 2 TeV. In Fig.2(l) there is a peak around V§ ~ 400 GeV . This peak comes 
from the ti threshold effects and enhanced by the top quark loop diagrams. 

The calculation shows the result of the cross section a is not sensitive to the masses of the 
related two neutral Higgs bosons m^o and m H o, especially when is far beyond the value 
regions of the two masses m h o and m#o. From the Feynman diagrams in Fig.l, we know that 
the diagrams, which are related to the neutral Higgs bosons and the trilinear Higgs couplings, 
are only appeared in the s-channel as shown in Fig.l(a,b), but their contribution is negligible 
quantitatively. This comes from two reasons. Firstly, the contribution is enhanced mainly 
by the Yukawa coupling which is proportional to m t /mw- Since the 7 exchange s-channel 
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diagram has no contribution and Z G boson exchange s-channel diagram does not include 
Yukawa coupling, the contribution from these kinds of s-channel diagrams is rather small. 
Among the s-channel diagrams, Fig. 1(b) has only one Yukawa coupling, whereas each box 
diagram in Fig.l(c ~ h) has two such couplings, therefore the cross section of the subprocess 
is mainly contributed by the box diagrams. Secondly, when the energy of CMS system \/I is 
not in the regions arround m^v, and m H o , the propagators in s-channel diagrams will suppress 
their contributions to the cross section. So it is clear that the measurement of this process 
could not be expected to give constraint on the features of neutral Higgs bosons and the 
trilinear Higgs couplings at the future pp colliders. 

In figure 3 we show the cross section a as the function of tan/3 with m,H± = 150 GeV, 
when \f§ = 400 GeV and V§ = 800 GeV respectively. The cross section goes down when 
tan/5 varies from 1.5 to 6, and raises up when tan/3 varies from 6 to 30. The cross sections 
are extremely small in the region arround tan (3 = 6. Fig. 4 is the plot of cross section a 
as the function of the charged Higgs boson mass with V§ = 1 TeV, tan/5 are 1.5 and 30 
respectively. There the cross section increases with increasing mass of the charged Higgs 
boson m,H±- In the region where m H ± < 300 GeV , the curves are flat, whereas they raise 
up quickly when m H ± > 300 GeV. That is because when the charged Higgs boson pair 
production threshold approaches the CMS energy \/J = 1 TeV, the cross section will be 
enhanced by the threshold effect of the charged Higgs boson pair production. 

In Fig. 5(1) and Fig. 5(2) we show the cross section of pp — > gg — > H + H~ + X process with 
m H ± being 150 GeV and 300 GeV respectively. The charged Higgs boson pair production 
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rates are read to be about 14 to 41 fermi-bar at the LHC energies when charged Higgs boson 
mass is 150 GeV and about 4.2 to 13.2 fermi-bar when m#± = 300 GeV. 

In conclusion, we have studied the pair production process of charged Higgs boson via 
gluon-gluon fusion in pp collider at LHC. The numerical analysis of its cross sections is 
carried out in the general two-Higgs- doublet model. With the possible parameters, the 
cross sections at future LHC collider can be 4 to 41 fermi-bar. These quantities are in 
the comparable range with that of the neutral Higgs boson pair production in gluon-gluon 
collisions as shown in Ref. |J. The calculation shows also that the production rate of charged 
Higgs boson pair is not sensitive to the neutral Higgs boson masses and the trilinear Higgs 
couplings, but it is strongly related to the couplings of the charged Higgs boson to fermions. 

This work was supported in part by the National Natural Science Foundation of China 
and the National Committee of Science and Technology of China. Part of this work was done 
when two of the authers, Ma Wen-Gan and Yu Zeng-Hui, visited the University Vienna under 
the exchange agreement (project number: IV.B.12). 

Appendix 



We adopt the same definitions of one- loop integral functions as in reference |TT]] and d 
4 — e. The integral functions are defined as 

(^r D f, D . i 



171 J 



[q 2 — m 2 ] 

{B -,B„;Br}(p,m 1 ,m i ) = —rj r -Jd g _ - - m2y 



[q 2 — m 2 ] [(q + p) 2 — m 2 ] 
{C ; C M ; C^; C fiup }(p 1 ,p 2 , mi, m 2 , m 3 
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( 2 W I d D 

in 2 J [q 2 - mf][(g + pi) 2 - m 2 ][(q + pi + P2) 2 - m 3 ] ' 

{.Do; -D^; A*"! D flupa }(p 1 ,p 2 ,P3, m 1 ,m2,m 3 ,m 4 ) = 

(27T/i) 4 D f ,D [}-'■! Qix'-iQiiv'-iQiivp'-iQuupa) 



S d " q 



in 2 J [q 2 - m 2 ][(q + p^ 2 - m 2 .]^ + px + p 2 ) 2 - ml][(q + px + p 2 + p 3 ) 2 - ml]' 
The form factors involved in equation (4) are represented in terms of the above integral 
functions explicitly as: 

/1 = 2m t (g h o H ± Ip g h o tt A h o + g H o H ± HT g H o tt A H o) d 5 [pi +P2, rra t , m t ] 

- 8m t (g h o H ± H ^g h o tt A h o + g H o H ± H ^g H o tt A H o)- (-4-1) 
(dC 24 + (p4 • P3)C )[-p 3 ,Pi + P2, "it, m t , m*] + ifcgl + ifi9+, 

where 

/f = fla + /f 6 + /fc, (-4 " 2) 



/f a = 2(2(m 6 m t C - m^C - m 2 H± C 12 ) + 2(pi • p 3 )(C + Cn - C12) - 

2(P2 -p 3 )(C + Cn) + 2(pi •p 2 )(Cn - C12) - dC 24 )[pi -p 3 , -p 4 ,m 6 ,m t ,m t ] 

+ 2(2(px • p 3 )(Cn - Co - C 12 ) - 2(px • p 2 )C n + 2(p 2 • p 3 )(C + 3C n - 2C 12 ) - 
dC 24 + 2m b (m t - m b )C - 4:m 2 H ±C n )[-p 2 + p 3 , -p 3 , m&, m t , m t ] 

+ A(2m 2 H ±(D 27 + 2D 3n - 2D 312 + 2L> 313 ) + 4m t (m t - m b )D 27 + 
2mf r± m t 2 (£>n - D 12 + D 13 ) + 2(pi • p 2 )(D 27 + 2£> 313 + m 2 D 13 ) + 
(pi • p 3 )(4 J D 312 - 4£> 313 - 2,D 27 + m 2 H± (D u - D 12 + 2D 21 + 2,D 22 + 
2D 23 - 4D 24 + 4L> 25 - 4D 26 ) + m b m t (2D 12 - D - 2D n - 2D 13 ) + 
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m 2 t (D + 2D n )) + 2(pi • p 2 ){pi • p 3 )(£> 23 + ^25 - A*) + 
2(p! • p 3 ) 2 (D 24 - D 22 - D 23 - D 25 + 2,026) + 
(P2 • p 3 )(m| ± (L>n - L> 12 + 2D 13 + 2D 23 + 2D 25 - 2D 26 ) - 
m b m t {D + 2D 13 ) + m 2 t {D + 2D 13 )) + 2(pi -p 2 )(p 2 • Pz){D 13 + D 23 ) + 
2(pi -P3)(P2 -Pz){D 2% - D 23 ))[-pi,pi -p 3 , -p4,m t ,m b ,m t ,mt\ 
+ A(Am 2 H± D 3U + 2m 4 H± (D u + D 21 + L> 24 ) - 4m b m t D 27 + 

4m t 2 D 27 - 2m b m t m 2 H ±(D + £> n + D 12 ) + 2m 2 H ±m 2 (D + 2£>n + L> i2 ) + 

2(pi • pa) (2D 312 + m 2 H± (D u + £> 12 + L> 21 + D 22 + 2D 24 ) - m b m t (D + D u + D 12 ) + 

m 2 (A) + D n + 2D 12 )) + 2( Pl ■ p 2 ) 2 (D 12 + D 22 + D 24 ) + ( Pl • p 3 )(2D 27 - AD 3l2 + 

4D 3 i3 + m 2 H± (2D 13 - D n - D 12 - 2D 21 - 2D 22 - AD 24 + AD 25 + 2D 26 ) + 

m b m t (D + 2D n + 2D 12 - 2L> 13 ) + m 2 (AD 13 - D - 2D n - AD 12 )) + 

2(pi • P2)(pi • Ps)(D 13 - D l2 - 2D 22 - 2D 24 + D 25 + 2D 26 ) + 

2(pi • p 3 ) 2 (£22 + ^23 + #24 - ^25 - 2D 26 ) + 

(pa • p 3 )(™|± (2D 25 - Dn - L>i2 - 4D 24 ) + m b m t (D + AD l2 - 2D 13 ) + 

m 2 (2D 13 -D - AD 12 )) + 2{ Pl ■ p 2 ){p 2 ■ p 3 ){D 2% - D 12 - 2D 22 ) + 

2(pi -P3)(P2 -P?){D 23 + 2£> 22 - 3D 26 ))[-pi, -p 2 +p 3 , -p 3 ,m t ,m b ,m t ,m t }, (A - 

fu = fia ( m b ++rnt), (A -A) 
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= 2d(B [-p 3 , m t , m t ] + B [-p 3 , m b , m b ] - B [-p 2 , m t , m b \) 

+ A(m 2 H± (C n + C 12 ) + ( Pl ■ p 2 )(C n + C 12 ) - ( Pl ■ p 3 )(2C + 3C n ) - 
(P2 -Pz){Cq + Cn + 2Ci 2 ))[-pi +P3, -p 2 ,m b ,m t ,m b ] 

+ 4{m 2 H± {C n -C - Ci 2 ) + (pi • p 2 )(C n - C - C12) - (pi • p 3 )Cn + 
(p 2 •P3)(2Ci2 - Co - 3Cii))[p 2 -P3, -p 2 ,m b ,m t ,m b ] 

+ A(2(m 2 b - 2m b m t + m 2 t )D 21 + 2m 2 H± m 2 (D 12 - D n ) + 2(pi • p 2 )(£>27 + m^D 
(pi • p 3 )(2D 27 + ml(D 11 - D 12 ) + m 2 H± (D 13 - 2D 21 - 2D 22 + AD 2A ) + 
m 2 (Do - D u - D 12 ) - m b m t D ) + 2(p x • p 2 )(pi • P 3 )(£>25 - D 2e ) + 
2(pi -pg) 2 ^ - ^24) + (P2 ■ p 3 )(m 2 H± (D 11 - D 12 + 2D 25 - 2^26) - 
m 2 b D 13 - m h m t Do + m 2 (D + D 13 )) + 2(p 1 ■ p 3 )(p 2 • p 3 )£> 2 6 - 
2(pi -p 2 )(p 2 -p 3 )(-Di 3 + £> 23 ))[pi, -pi +p 3 , -p 2 ,m t ,m b ,m u m b ] 

+ 4(2(m^ + m^± + m 2 )D 27 - Am b m t D 27 + m 2 b m 2 H± (D u + D 12 - D 13 ) + 
m%± (D 13 - D n - D 12 - 2D 21 - 2D 24 + 2D 25 ) + 
m 2 H± m 2 {D 13 - 3D n - D 12 ) + ( Pl ■ p 2 ){4D 27 + m 2 b {D u + D 12 - D 13 ) + 
m 2 H± (D 13 - D u - D 12 - 2D 21 - 2D 22 - 2D 23 - AD 24 + AD 25 + AD 2% ) + 
m 2 (3D 13 - D u - 3D 12 )) + 2{ Pl ■ p 2 ) 2 (D 25 - D 22 - D 23 - D 2A + 2D 2% ) + 
(Pi • Ps)K 2 (^o + D vl + 3Di 2 ) - m b m t D - m 2 b (D n + D 12 ) + 
m 2 H± (2D U + 2D 12 - D 13 + 2D 21 + 2D 22 + AD 2A - 2D 26 - 2D 27 )) + 
2(pi • p 2 )(pi • p 3 )Pi 2 - D 13 + 2D 22 + 2L> 2 4 - ^25 - 2D 26 ) - 
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2(pi • p 3 ) 2 (D 12 + 22 + ^4) + (P2 • P3)(mg(£>i3 - 2Di 2 ) - m 6 m t L> + 
m 2 H± (D n + 0> 12 + 40 24 - 20 25 ) + m t 2 (O + 20 12 - D 13 )) + 
2(pi • p 2 )(p 2 • p 3 )p23 + 20 22 - 30 26 ) + 

2(pi ■ Ps)(P2 ■ P3)(D 26 - D 12 - 2D 22 ))[p 1 ,p 2 -p 3 , -p 2 ,m t ,m b ,m t ,m b ], (A - 5) 



/+ = /fK^-m 6 ]. (A -6) 



/ 2 = Sm t (g h o H ± H ^g h o tt A h o + g H o H ± H ^g H o tt A H o)- (A — 7) 

(C u + 2C 21 + 4C 22 - 4:C 23 )[-p 3 , Pl + p 2 , m t , m t , m t ] + if 2 g 2 _ + if 2 + g 2 + , 1 > 

where 

f2 = /£ + & + / 2 "c, {A " 8) 

f2a = 4 (C + C u - C12) [pi - P3, -P4, m 6 , m t , m t ] 
+ 4(C + Cn - C12) [-p 2 + P3, -P3, m 6 , m t , m t ] 

+ (8(0 27 + 4L>3ii - 2D 3 i2 + 2D313) + 4m^ ± (0> 12 - D u - AD 21 - 2D 22 + QD 24 - 
40 25 + 20 26 - 40 3i + 8O310 + 8O34 - 8O35 - 40 36 - 40 37 ) + 

Am b m t {D + 40> n - 40 12 + 40 13 + 40 2i - 40 24 + 40 25 ) + 4m 2 (40 12 - D - 40> n - 
40 13 - 40 2i + 40 24 - 40 25 ) + 8(pi • p 2 )(D 22 - 24 - 26 + 2D 3W - 2D 35 - 20 37 ) + 
8(pi • p 3 )(0 2 2 + 23 - 24 + 25 " 2,026 - 4O 3 i - 20 34 + 20 35 + 20 36 + 20 37 ) + 
8(^2 • Pa) (^24 - 22 - 23 - 25 + 20 26 )) [-pi, pi - p 3 , -P4, m t , m b , m t , m t ] 
+ 4(2(0 27 + 4O311 - 2O312 + 2O313) + m 2 H± (0 12 - O n - 40 2i - 20 22 + 

60 24 - 40> 25 + 20) 26 - 4D 3 i + 40 34 - 40> 35 ) - m t (m t - m b )(D + 40 n - 40 12 + 
40 13 + 40 2 i - 40 24 + 40 25 ) + 2(pi • p 2 ){D 22 - 24 - 26 - 2O 3 i - 20 34 + 20 36 ) + 
2(pi • p 3 )(0 24 - 22 - 23 - 25 + 20 26 + 4O310 + 2O34 - 2O35 - 20 36 - 20 37 ) + 
2(P2 -P3)(0 2 2 + 23 - 24 + 25 - 20 26 ))[-pi, -P2 +P3, -p 3 ,m t ,m b ,m t ,m t ], 

(A -9) 



f2b = f2a (m b ^m t ), 04-10) 
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Ac = -4Cn[-pi + p 3 , -p 2 , m b , m t , m b ] - AC U [P2 - p 3 , -p 2 , m b , m t , m b ] 

+ A{m 2 b (D 12 - D u - 2D 21 + 2D 24 ) + m 2 (2D 24 - D - 3D U + 3D 12 - 2D 21 ) + 
m b m t (D + AD U - AD 12 + 4^21 - AD 2A ) + m 2 H± (2D 26 - D 13 - 2D 25 ) - 2D 27 + 
2(pi • p 2 )(D 22 - D 24 ) + 2(pi • p 3 )(£> 2 4 - £22) + 
2(P2 -P3)(^24 - #22 - Ae))[Pi, -Pi -p 2 , m t , m b , m t , m b ] 
+ A(m 2 b (D 12 - D u - 2D 21 + 2D 24 ) - 2D 27 + m 2 H± (2D l2 - D 13 + AD 24 - 2D 25 ) + 

m b m t (D + AD U - AD 12 + AD 21 - AD 24 ) + m 2 (3D 12 - D - 3D U - 2D 21 + 2D 2A ) + 
2(pi • p 2 )( J Di 2 + £22 + D 2A - D 26 ) - 2(pi • p 3 )( J Di2 + ^22 + + 
2(p 2 -p 3 )(-D 2 6 - -D12 - -D22 - ^24))[pi,P2 -p 3 , -p 2 , m t , m b , m t , m b ], 

(A -11) 

f+ = f~[m h ^-m b ]. {A -12) 



f 3 = Sm t (g h o H ± H ^g h o tt A h o + g H o H ± H ^g H o tt A H o)- (A -13) 

(C u + 2C 21 + AC 22 -AC 23 )[-p 3 ,p 1 +p 2 ,m t ,m t ,m t }+if 3 g 2 _ + if 3 + g 2 + , 1 ' 

where 

/»" = & + £ + (i* - 14) 

/3a = A(C - C 12 )[-p 2 +p 3 ,-p 3 ,m b ,m t ,m t ] 

+ 4(2C 12 - 3C - AC U - 2C 21 + 2C 23 ) [pi - p 3 , -p 4 , m b , m t , m t ] 
+ A(A(D 27 + D 312 ) + m 2 H± (D 12 - D n - AD 21 - 2D 22 - AD 23 + 6D 24 - 
8D 25 + 6^26 - 8^3X0 - 4L> 32 - AD 34 + 8D 36 + 8^ - AD 39 ) + 
m b m t {D + ADu - AD l2 + AD 13 - AD 22 + AD 24 + AD 26 ) - m 2 {D + AD U - 
AD 12 + 4D 13 - 4£> 22 + 4L> 24 + 4£> 26 ) + 2( Pl ■ p 2 )(D 22 - 2D 23 - D 24 - 2D 25 + 
D 26 - 2D 310 + 2D 38 - 2D 39 ) + A( Pl ■ p 3 )(D 22 + D 23 - D 24 + D 25 - 
2D 2% + D 310 + £> 32 - D 3& - 2£>38 + D 39 ))[-p 1 ,p 1 -p 3 , -p4,m t ,m b ,m t ,mt] 
+ A(2(2D 313 - D 27 ) + m 2 H± (2D 26 - D n - D 12 - 4L> 25 - 4L> 35 )- 

m t (m t - m b )(D + 4D 13 + AD 25 ) - 2( Pl ■ p 2 )(D 12 + D 26 + 2D 310 ) + 

2(pi • Ps)(D 12 - D 13 - D 23 + D 26 + 2D 310 - 2D 37 ) + 

2(P2 -P 3 )(-Di2 - D 13 + £> 23 - D 26 ))[—pi, -p 2 +p 3 , -p 3 ,m t ,m b ,m t ,m t }, 

(A -15) 



/ 3 6 = / 3a (m b ^m t ), (A - 16) 
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fsc = 4(C - Ci 2 ) [P2 - P3, -P2, m b , m t , m b ] + A(C 12 + 2C ) \-p l + p 3 , -p 2 , m b , m t , m b ] 
+ A(m 2 (D 12 - D n + 2D 22 - 2D 24 + 2D 25 - 2D 26 ) + m 2 H± (2D 21 - D 13 + 
2D 22 - AD 2i ) - 4D 27 + m b m t (D - AD 22 + 4D 24 - 4L>2 5 + 4D 2 6) + 
m 2 (D u -D - D 12 + 2D 22 - 2D 24 + 2D 25 - 2D 26 ) + 
2(pi -P2)(D 22 - D 24 ) + 2( Pl -p 3 )(D 26 - D 25 + 2D 24 - 2D 22 ) + 
2(P2 ■P3)(D 13 + D 23 - D 26 ))[p u -p x +p 3 , -p 2 ,m t ,m b ,m t ,m b ] 
+ A(m 2 H± (2D 13 + 2D 25 - D u - D 12 ) - 2D 27 + m 2 b (D 13 + 2D 25 ) + 

m b m t (D - AD l3 - AD 2b ) + m 2 (2D 25 - D + 3£> 13 ) + 2( Pl ■ p 2 )(D 25 + D 26 - 

D 12 + 2D 13 - D 23 ) + 2(pi • p 3 )(D 12 - D 13 - D 25 - D 26 ) + 

2(P2 -P3)(D 12 - D 13 + D 23 - D 26 ))\p!,p 2 -p 3 , -p 2 ,m t ,m b ,m t ,m b ], 

(A - 17) 

f^f.lm.^-m,]. (A -18) 



f 4 = 8mt(gh°H±H^gh"ttA h o + gH°H±H*gH°ttAH°)- 

(C11 + 2C 2 i + 4C 22 - 4C23) [-P3, Pi + P2, mt, m, m t \ + if 4 9- + if£g\, 



(A - 19) 



where 



J! = fla + fib + flc, {A - 20) 



ha = 4 (C - C12) [pi - p 3 , -Pa, mb, m t , m t \ 

+ 4(2C 12 - 3C - 4C n - 2C 2 i + 2C 23 )[-p 2 + P3, -P3, m b , m t , m t ] 
+ A(2(2D 313 - D 27 ) + m 2 H± (D 12 - D n - 2D 13 - AD 25 + 2D 26 + 

AD 310 - AD 35 - AD 37 ) + m b m t (D + AD 13 + AD 25 ) - m 2 (D + AD 13 + AD 25 )- 
2(pi • P2)(D 13 + D 26 + 2D 37 ) + 2(pi • p 3 )(D 23 - D 26 - 2D 310 + 2L> 37 )+ 
2(P2-P3)(D 26 - D 23 ))[-pi,pi - p 3 , -p 4 ,m t ,m b ,m t ,mt\ 
+ A(A(D 27 + D 312 ) + m 2 H± (3D 12 - D u - 2D 13 - AD 21 + 2D22 + 2D 2 4- 

4£>2 5 - 2^26 - 4^310 - 4£> 3 4 + 4^36) - m t (m t - m b )(D + AD U - AD 12 + 
AD l3 - AD 22 + AD 24 + AD W ) + 2{ Pl ■ p 2 )(D 12 - D l3 + 3D 22 - 3D 24 - 3D 26 + 
2D 32 - 2D 36 - 2D 38 ) + 2(p! • p 3 )(D 13 - D 12 - 2D 22 - 2D 23 + 2^24 - 2^25+ 
4^26 - 2^310 - 2D 32 + 2^36 + 4,038 - 2D 39 ) + 
2(P2 -P3)(D 13 - D 12 ))[-p u -p 2 +p 3 , -p 3 ,m t ,m b ,m t ,m t }, 

(A -21) 



Ub = f4a (m b ^m t ), (A -22) 
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/ 4c = 4(C - C 12 )[-pi + p 3 ,-p 2 ,m b ,m t ,m b } + A(C 12 + 2C ) [p 2 - Ps, ~P2, m&, m t , m b ] 
+ A(m 2 H± (D 12 - D u + 2D 23 ) - 2D 27 + m 2 (D 13 + 2D 25 ) + 

m b m t (D - AD 13 - AD 25 ) + m 2 t (2D 25 -D + 3D 13 ) + 2( Pl ■ p 2 )(D 13 + D 26 ) + 

2(Pi ■P3)(D 25 - D 26 ) - 2(p 2 -p 3 )(D 23 + D 26 ))\p u -pi+p 3 , -p 2 ,m t ,m b ,m t ,m b } 
+ A(m 2 (D 12 - D u + 2D 22 - 2D 24 + 2D 25 - 2D 26 ) + 

m 2 H± (2D 21 - D l3 + 2D 22 + 2D 23 - 2D 25 - AD 26 )+ 

m b m t (D - AD 22 + 4D 24 - AD 25 + AD 26 ) - AD 27 + 

m\{D lx -D - D 12 + 2D 22 - 2D 24 + 2D 25 - 2D 26 ) + 

2(pi • p 2 )(D 22 + D 23 + D 24 - D 25 - 2D 26 ) + 2( Pl ■ p 3 ){D 13 - 2D 24 + D 25 + D 26 )+ 
2(p 2 - p 3 )(SD 26 - 2D 22 - D 23 ))[p 1 ,p 2 - p 3 , -p 2 ,m t ,m b ,mt,m b ], 

(A - 23) 

ft = fr[m b ^-m b ]. (A -24) 

f 5 = 8m t (ghOH±H? i gh t HtAhO + gHOH±H^9HOttAHo)- (A_OK\ 

(C n + 2C 21 + AC 22 - AC 23 ) [-p 3 , Pl + p 2 , m t , m t , m t ] + %^ g 2 _ + tf+g 2 + , {/i Z0) 

where 

ft = fsa + & + &, (A - 26) 

fba = 4 (Co + C u + 2C 12 + 2C 23 ) [pi - p 3 , -p 4 , m b , m u m t \ 
+ A(C + Cn + 2C 12 + 2C23) [-P2 + Ps, -P3, m b, m t , m t ] 

+ (Am 2 H± (D 12 - D n - 2D 13 - AD 23 - AD 25 + 2D 26 - AD 310 + AD 38 - AD 39 )- 
Am t (m t - m b )(D + AD l3 + AD 26 ) - 8( Pl ■ p 2 )(D 13 + 2D 23 + D 26 + 2D 39 ) + 
16(pi -P3)(D 23 - D 26 - D 38 + D 39 ))[-p 1 ,p 1 -p 3 , -p 4 ,m t ,m b ,m t ,m t ] 

+ A(m 2 H± (D 12 - D u - 2D 13 - AD 25 - 2D 26 - AD 3W )- 

m t (m t - m b )(D + AD 13 + 4^20) - 2( Pl ■ p 2 )(D 13 + 3D 26 + 2D 38 ) + 
4(pi -Pz){D 2 % - D 23 + D 38 - ^ 39 ))[-pi, -p 2 +p 3 , -p 3 ,mt,m b ,m t ,mt\, 

(A - 27) 



hb = ha (m b ^m t ), (A -28) 
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= 4(C + Cu) [-pi + p 3 , -p 2 , m b , m t , m b ] + 4(C + C n )[p 2 - p 3 , ~P2, m b , m t , m b ] 
+ A(m 2 (D 13 - 2D 23 + 2D 26 ) + m 2 H± (D 12 - D u - 2D 25 + 2D 26 ) + 

m b m t {D + AD 23 - AD 26 ) + m 2 t {2D 26 - D - D 13 - 2D 23 ) + 

2(pi -p 2 ){D 13 + D 26 ) - 2(pi -p 3 )(D 13 + L>26))[pi, -Pi +P3, -p 2 ,m t ,m b ,m t ,m b } 
+ 4(m 2 H± (D 12 - D n - 2D 13 - 2D 25 ) + m 2 (D 13 - 2D 23 + 2D 26 ) + 

m b m t (D + AD 23 - AD 26 ) + m 2 (2D 26 - D - D 13 - 2D 23 ) + 

2(Pi -Pz){D 13 + D 26 ))\p u p 2 -p 3 , -p 2 ,m t ,m b ,m t ,m b }, 

(A - 29) 

U = f 5 '[m b ^-m b ]. (A - 30) 



h = ha + hb + he, (A - 31) 

= 9-9+(8(C 12 - Cn)[pi -p 3 , -p 4 ,m t ,m b ,m b ] - 8C u [-p 2 + p 3 , -p 3 ,m b ,m t ,m t ] 
+ 16((pi • p 3 )(D 12 - D n - D 13 - D 21 - D 22 - D 23 + 2D 24 - 2D 25 + 2D 26 ) + 

{P2-p-i){D 2& - D 13 - D 23 - D 25 ) - D 27 )[-p u pi -p 3 , -p 4 ,m t ,m b ,m t ,mt\ 
+ lQ(m 2 H± (D 12 + D 22 + D 24 ) - 2D 27 + ( Pl ■ p 2 )(D 12 + D 22 + D 24 ) + 

(pi • p 3 )(D 26 - Dn - D 21 - D 22 ) + 

(P2-P3)(D 26 - D 12 - D 22 - D 24 ))[-p u -p 2 +p 3 , -p 3 ,m u m b ,m u m t \), 

(A - 32) 

hb = ha (m b ^m t ), (A - 33) 



= ^9s9zH±H^9ztt A z ■ {dB^-pi - p 2 , m t , m t \ - AC 24 [p 4 , -p x - p 2 , m t , m t , m t ]}+ 

g-g + (8(C 12 - Cn)[-pi + p 3 , -p 2 ,m b ,m t ,m b ] 
+ 8C + Cu - C 12 )\p 2 - p 3 , -p 2 , m b , m t , m b ] 

- \&D 27 \pn -pi +p 3 , -p 2 ,m t ,m b ,m t ,m b ] + 8(m 2 H± (D n - D 12 + D 13 + 2D 21 - 
2D 24 + 2D 25 ) + m 2 (D 12 - D n - D 13 ) + m 2 (D n - D 12 + D 13 ) - AD 27 + 
2(pi • p 2 )(D 24 - D 22 - D 23 - D 25 + 2,026) + 
2(Pi -Pz){D 22 - D 24 - D 26 ))[p u p 2 -p 3 , -p 2 ,m u m b ,m u m b \). 

(A - 34) 



h = ha + hb + he (A - 35) 
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ha = g-g+(S(C + C n - C 12 )\pi - p 3 , -Pa, m b , ™t, m] 
+ 8(C + C n - C 12 ) \-p 2 + P3, ~P3, m b , m t , m t ] 

+ 8(m 2 H± (D n - D 12 + 2D 13 + 2D 21 + 2D 22 + 4L> 23 - AD 24 + QD 25 - QD 26 )- 

2D 27 - m 2 t D + 2(pi • p 2 )(D u - D 12 + D 13 + D 21 + D 22 + 2D 23 - 2D 24 + 3^25- 
3,026) + 2(p 4 ■ p 3 ) (L> 2 4 - ^22 - D 23 - D 25 + 2D 26 )) [-p u p 1 - p 3 , -p 4 , m t , m b , m t , m t ] 

+ 8{m 2 H± (D u - D 12 + 2D 13 + 2D 2 i - 2D 22 + 2D 25 + 2D 26 ) - 2D 27 - m 2 t D + 
2(pi • p 2 )(D 11 - D 12 + D 13 + D 21 - D 22 + D 25 + £ 26 ) + 

2{Pi-p 3 ){D 22 + £>23 - £>24 + 0> 25 - 2L>2 6 ))[-Pl, ~P2 + P3, ~P3, "l t , 77l 6 , m t , m t ]) . 

(A - 36) 



hb = ha (m b <-*mt), (A - 37) 



/ 7c = 2ig 2 g ZH ± H ^g ztt A z ■ {-dB [-pi - p 2 , m t , m t ] + 

A[(m 2 H± + pi ■ p 2 - pi ■ p 3 - p 2 ■ p 3 )Cn + dC 2 4][P3, ~Pi - P2, "it, "it, mt]} 
+ g-g+(-%( C n + 2C )[-pi+p 3 , -p 2 ,m b ,m t ,m b ]- 

8(Cn + 2C ) [P2 - P3, -P2, m b , m t , m b ] 
+ 8{m 2 H± (D 13 - 2D 2 i - 2D 22 + 4L> 24 ) + m 2 (D u - D 12 ) + QD 27 + 

m 2 {D - D u + D 12 ) + 2(p! • p 2 )(L> 25 - £> 26 ) + 2(pi • p 3 )(£> 2 2 - D 24 )- 

2(p 2 • Pz)Di 3 ) [pi, -pi + p 3 , -p 2 , m*, m b , m t , m b ] 
+ 8{m 2 H± (2D 12 - D 13 - 2D 21 + 2D 24 ) + m 2 (D - D u + D 12 ) + m 2 (D u - D 12 )+ 

QD 27 + 2(pi • p 2 )(L»i 2 - D 13 + D 22 - D 24 + D 25 - D 26 ) + 2( Pl ■ p 3 )(D 24 - 

D 12 - D 22 ) + 2(p 2 -p 3 )(D 13 - D 12 ))[p 1 ,p 2 -p 3 , -p 2 ,m t ,m b ,m t ,m b ]). 

(A - 38) 



/8 — /8a + hb + he, 

ha = g-g+(8(C - Cn)[pi -p 3 , -p 4 ,m b ,m u m t ] 
+ 8(C - Cn)[-p 2 +p 3 , -p 3 ,m b ,m t ,m t ] 

+ 8(m 2 H± (D 12 - D n + AD 23 + 2D 25 - 2D 26 ) - 8D 27 - m 2 D + 

2(pi -p 2 )(D 25 + 2D 23 - D 26 ) + 2( Pl -p 3 )(D 24 - D 22 - 2D 23 - D 25 + 
3D 26 ) +4(p 2 -P3)(-D 2 6 - D 23 ))[-pi,pi - p 3 , -p 4 ,m t ,m b ,m t ,mt\ 

+ 8{m 2 H± (D 12 - D u + 2£> 2 2 + 2D 24 ) - 8D 27 - m 2 D + 

2(pi -P2)(^22 + £24) + 2(pi -p 3 ){D 2b + D 26 - D 22 - D 24 ) + 
4(P2 •P3)(0 ) 26 - D 22 ))[-pn -p 2 +p 3 , -p 3 ,m t ,m b ,m t ,m t ]), 



(A - 39) 



(A - 40) 
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(A - 41) 



he = 2iglg ZH ± H ^g a ztt A z ■ {dB [-pi - p 2 , m t , m t } + 

4[(pi -P3+P2-P3-Pi-P2- m 2 H± )C n - dC 2A ][p 3 , ~Pi - P2, m t , m t , m t }} 
+ g-g+(-8(C + Cn)[-pi + p 3 , -p 2 ,m b ,m t ,m b ] 
- S(C + C u )\p 2 -P3,-P2,m b ,mt,m b ] 

+ 8(m 2 H ±(D u - D 12 + 2L> 25 - 2D 26 ) + m 2 t (D + D 13 ) - m 2 D l3 - 

2(pi -p 2 ){D n + D 23 ) + 2(pi -p 3 )(-Di3 + £>2e))[pi, -Pi +P3, -p 2 ,m t ,m b ,m t ,m b ] 

+ 8(m^ ± (L> n - D 12 + 2D 13 + 2D 25 ) + m?(A) + A3) - ™^i3+ 

2(pi -P2)(D 26 - D 23 ) - 2(pi -p 3 )(-Di3 + Ae))[pi,P2 -P3, -p 2 , m t , m 6 , m t , m 6 ]). 

(A - 42) 



/ 9 = /9a + /% + /9c, (A - 43) 

/ 9a = 16^_^ + ((Di2 - D n - D13 - D 21 + L> 24 - D 25 )[-p 1 ,p 1 -p 3 , -p 4 ,m t ,m b ,m t ,m t ] 
+ (D 12 - D u - D 13 - D 21 + D 24 - D 25 )[-pi, -p2 +p 3 , -p 3 ,m t ,m b ,m t ,m t ]), 

(A - 44) 



fab = ha (m b <r+m t ), (A - 45) 

f 9c = 32ig 2 s g ZH ± H Tg Ztt Az ■ (C 23 - C 22 ) [p 3 , -p x - p 2 , m t , m t , m t \ 

+ 16g-g+((D 22 - D 24 + D 25 - D 26 )\p 1 ,-p 1 +p 3 ,-p 2 ,m t ,m b ,m t ,m b \ (A - 46) 
+ (D 13 + D 25 )\pi,p 2 -p 3 , -p 2 ,m t ,m b ,m t ,m b ]). 

ho = flOa + /106 + hoc, (A ~ 47) 

hoa = 16g-g+({D u - D 12 + D 13 + D 21 - D 24 + D 25 )[-p 1 ,p 1 - p 3 , -p 4 ,m t ,m b ,m t ,m t ] 
+ (D n - D 12 + D 13 + D 21 - D M + D 25 )[—pi, -p 2 +p 3 , -p 3 ,m t ,m b ,m t ,mt\), 

(A - 48) 



hob = hoa (m b ^m t ), (A - 49) 
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/ioc = lQig 2 s g Z H±H^9ztt A z • (Cii + C21 + 2C 22 - 2C 23 )[p 4 , -pi - P2, m t , m t , m t ] 
+ 16^_g- + (-(Di 3 + D 25 )\pi, -pi +p 3 , -P2,m u m b ,m u m b ] 
+ (L>24 - ^22 - -D25 + D 2 q)[pi,P2 -P3, -P2,m t ,m b ,mt,m b ]). 

(A 



fll — flla + fllb + flic, (A 

f lla = 16g^g + ((D 2 2- D 13 - D 2 4- D 2 6)[-pi,-p2+P3,-P3,mt,m b ,mt,mt} 



(D 13 + D 25 )[-pi,Pi -P3, -P4,m t ,m b ,m t ,mt\), 



fua = 16g-g+((D 13 - D 22 + D 2i + D 26 )[-p 1 ,p 1 - p 3 , -p 4 ,m t ,m b ,m t ,m t } 
+ (D 13 + D 25 )[-p u -p 2 +p 3 , -p 3 ,m t ,m b ,m t ,m t ]), 



(A 



fiu = fiia (m b <-> m t ), (A 

f Uc = 32ig 2 s g ZH±H ^g a Ztt A z ■ (C 23 - C 22 )[p 3 , -p x - p 2 , m u m t , m t ] 

+ l6g-g+((D 2 Q - £> 23 )[pi, ~P\ +P3, -P2, m t ,m b ,m t ,m b ] (A 
+ (D 26 - D 23 )\p 1 ,p 2 -p 3 , -p 2 ,m t ,m b ,m t ,m b ]). 



fl2 — fl2a + fl2b + fl2c, (A 

(A 



fi2b = f\2a (m b <-> m t ), (A 

/12c = 16ig 2 gzH±H^g a zttAz ■ (C u + C 2 i + 2C 22 - 2C 23 )[P4, -Pi - P2, m t , m t , m t ] 
+ ^6g-g+((D 23 - £> 26 )[pi, -Pi +P3, -p 2 ,m t ,m b ,m t ,m b ] 
+ (D 23 - £> 26 )[pi,p 2 -p 3 , -p 2 ,m t ,m b ,m t ,m b ]). 

(A 

In the above form factor expressions we used the following definitions. 

-2i 
9jtt = -g-e, 

yztt 4cos# v 3 ; ' 
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a W 

9ztt ~ 4cos0' 



g^HT = -ie, 

—ig cos 29 



9ZH ± Hf — ■ 

9HHt = -ig 
9h°tt = -ig 



2cos# ' 
m t sin a 

2mw sin (3 ' 
m 4 cos a 

2mw sin /3 ' 



• / In \ cos 2/3 

9h°h±h* = -ig{m w cos{P - a) - m z - — cos(/^ + a)), 

2 cos 9 W 

• / • In \ COS 2/3 . . _ 

gh?H±H^ = -ig{m w sm((J -a)+ m z - — sin(/3 + a)), 

Z COS (7w 

g + = — y= {m b tan/3 + m t cot (3), 

ig 

g_ = — -= (m^ tan f3 — m t cot (3), 

2\/2mw 

A - — 

s ' 

. — i 

^ = 2"' 

s — 

= 



s — m 2 H0 + iT H om H o ' 

A„o = r -^- . (A - 59) 

s - m z h0 + iT h om h o 

The arguments of two-point, three-point and four-point integral functions are written at 
the end of formulae in paranthesis. The repalcement m b — > — m fe is not performed on the 
arguments of loop intergral functions B, C and D. 
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Figure captions 

Fig.l The Feynman diagram of subprocess gg — > H + H~ . (a,b,c,d)s-channel diagrams, (e ~ 
j)box diagrams including t- and u-channels. 

Fig.2(l)Total cross sections of the subprocess gg — > H + H~ as the functions of the with 
m h o = 100 GeV, m H a = 150 GeV and m H ± = 150 GeV. The full line is for tan/5 = 1.5. 
The dashed line is for tan (3 = 30. 

Fig.2(2)Total cross sections of the subprocess gg — > H + H~ as the functions of the \/§ with 
m h o = 100 GeV, m H o = 150 GeV and m H ± = 300 GeV. The full line is for tan/3 = 1.5. 
The dashed line is for tan j3 = 30. 

Fig.3Total cross sections of the subprocess gg — > H + H~ as the functions of the tan/3 with 
m h o = 100 GeV, m H o = 150 GeV and m H ± = 150 GeV. The full line and dashed line 
correspond to \/J = 400 GeV and = 800 GeV respectively. 

Fig.4Total cross sections of the subprocess gg — > H + H~ as the functions of the m H ± with 
m h o = 100 GeV, m H o = 150 GeV and y/I = 1 TeV. The full line is for tan/3 = 1.5. The 
dashed line is for tan/3 = 30. 

Fig.5(l)Total cross sections of the process pp — > — > H + H~ + X as the functions of 
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the y/s with m h o = 100 GeV, m#o = 150 GeV and m H ± = 150 GeV. The full line is for 
tan (3 = 1.5. The dashed line is for tan/5 = 30. 

Fig.5(2)Total cross sections of the process pp — > gg — > H + H~ + X as the functions of 
the v/i with m^o = 100 GeV, m^o = 150 GeV and m H ± = 300 GeV . The full line is for 
tan/3 = 1.5. The dashed line is for tan/5 = 30. 
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